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Abstract 

The Morton- or z-curve is one example for a space filling curve: Given a 
level of refinement L S No, it maps the interval [0, 2‘^^) DZ one-to-one to a 
set of J-dimensional cubes of edge length 2“^ that form a subdivision of the 
unit cube. Similar curves have been proposed for triangular and tetrahedral 
unit domains. In contrast to the Hilbert curve that is continuous, the 
Morton-type curves produce jumps. 

We prove that any contiguous subinterval of the curve divides the do¬ 
main into a bounded number of face-connected subdomains. For the hy¬ 
percube case in arbitrary dimension, the subdomains are star-shaped and 
the bound is indeed two. For the simplicial case in dimensions 2 and 3, 
the bound is proportional to the depth of refinement L. We supplement 
the paper with theoretical and computational studies on the frequency of 
jumps for a quantitative assessment. 


1 Introduction 

The Peano curve [17] and the Hilbert curve [13] are continuous maps from the 
line onto the d-dimensional unit cube. A large number of such space filling curves 
(SFC) has been described in the literature; see for example [3,12, 18] and the 
references therein. They are usually defined in terms of a recursive prescription. 
For numerical applications, the curve is made discrete and finite by bounding 
the depth of the recursion. The smallest units of space that are traversed may be 
called elements. The Morton- or z-curve, originally described by Lebesgue [14] 
and adapted to data storage by Morton [16], also creates such a map, but it is not 
continuous. In fact, it contains jumps throughout its length (see Figure 1). This 
raises the concern that a subsection of the curve may divide the space covered 
by its image into a large number of disconnected subdomains. Especially when 
the curve is used to divide a computational mesh between different processors 
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Figure 1: Contiguous subsections of the Morton curve at refinement level L = 3. Observe 
the jumps when the z-curve runs diagonally. In the right hand image, this produces two 
disconnected subdomains. In both pictures shown, the domain decomposes into two 
star-shaped pieces. 


for parallel computation, see e.g. [1,2,6,11,20], the surface-to-volume ratio of a 
fragmented subdomain could grow without bounds, which would likely increase 
the amount of data to be communicated. 

In this paper we eliminate that concern by proving that the (classic cubical) 
Morton curve can lead to no more than two subdomains, where we define a set of 
elements to be of the same subdomain if they are connected by a finite number 
of element face connections: 

Theorem 1. A contiguous segment of a Morton curve through a uniform or 
adaptive tree of maximum refinement level L produces at most two distinct face- 
connected subdomains. This result is independent of the space dimension. 

We note that a proof for the two-dimensional square has been given in [3, 
pages 175-177] that proceeds by illustrating and enumerating a finite number of 
cases. (In fact, we adapt these ideas to dimensions two and three in Section 3, 
and also restate the extension to adaptive meshes in Section 4.3.) It is said in [3] 
that the construction extends to dimensions three and higher. This is entirely 
plausible, yet we see that the number of cases to discuss grows with the space 
dimension and would eventually require some kind of automation. Thus, we 
proceed inductively over d to provide dimension-independent results. We also 
supply a formal non-inductive proof to show that the connected segments are 
star-shaped. 

For the triangular and tetrahedral Morton curves introduced recently [7], we 
show that the bound is proportional to the depth of refinement L: 

Theorem 2. A eontiguous segment of a tetrahedral Morton eurve through a 
uniform or adaptive tree of maximum refinement level L > 2 produces at most 
2{L — 1) faee-connected subdomains in 2D and at most 2L-\-l in 3D. For L = 1 
there are at most two face-eonnected subdomains. 
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We complete our study with a statement on the lower bound on the fraction of 
continuous segments and provide an algorithm and numerical results to illustrate 
the distribution of continuous vs. discontinuous segments. This supports our 
conjecture that the tetrahedral Morton curve is no worse in practice than the 
original cubical construction. 

2 Concepts and notation 

There is a natural identification between Morton-ordered elements on the one 
hand and uniform and adaptive quadtrees [10] and octrees [15] on the other. 
This is true for the tetrahedral Morton curve [7] as well. We will often refer to 
the elements as (sub)quadrants irrespective of the shape or space dimension d. 
Different ways exist to formalize the definition of a general space filling curve; 
one is to identify a finite set of types of transformations and rules to apply them 
recursively [12]. In this document we restrict the theory and notation to the 
minimum required to treat the cubical and the tetrahedral Morton curve. 

2.1 The cubical Morton curve 

The Morton subdivision of a d-dimensional hypercube [16] can be constructed 
by recursion. When dividing a cube into 2^^ half-size subcubes, we enumerate 
these with the binary index 

g= (grf...gi)2 G [0,2'^)nZ (1) 

comprised of d bits qi G {0,1}. (We will drop PlZ in the following when it is clear 
that we are referring to integers.) Each of the bits i corresponds to the position 
of that subcube in the Xi coordinate direction, where 0 denotes the lower and 
1 the higher half. In our convention the most significant bit corresponds to the 
last dimension (z in three dimensions) and the least significant bit to x = xi. 
When counting through the possible values of q we see that the x\ coordinate 
changes its value fastest and the Xd coordinate slowest. Before a bit at position 
i flips, all numbers in in the lower i — 1 bits have to be counted through first. 

We can state one central and well known fact at this point: The flip of the ith 
bit amounts to a shift of the corresponding subcube parallel to the coordinate 
direction i. If we flip from zero to one, we move up, and else we move down 
the axis. It is easy to see that flipping one bit transforms the subcube into its 
neighbor across a face with normal direction ±Xj. 

We define a recursion by subdividing each subcube further using the same 
prescription. The root cube is associated with level £ = 0, with levels increasing 
with each subdivision. Subcubes exist at any level £ and are identified with the 
root of a corresponding subtree. Level-L subtrees are also called subquadrants. 
We count the sequence of level L (sub) quadrants with the index 

Q= (gl...g^)2 G [0,2'^^), (2) 

where each level-wise index q^ is defined as in (1). They designate the choice of 
subquadrants from the first subdivision = 1 to the last at level £ = L. This 
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sequence of choices can be understood as the path from the root to the leaf of a 
decision tree, where each decision is between 2'^ possibilities. The subset of 
occupied by the quadrant with index Q is 

Q(Q) ■.=[2-^{q\ql... gf ) 2 , ... gf )2 + 1)] x 

[2~^{q\ql... g 2 ) 2 , 2"-^((g^g|... g ^)2 + 1)] x 

: 

■ ■ ■ (Idh, ‘^~^{{QdQd ■ ■ ■ Q-dh + !)]■ 

We define a full or complete subtree by the set of all its descendant quadrants. 
A subtree is incomplete if the quadrants form a strict subset of descendants that 
are contiguous with respect to the indexing (2). We call such a subset a segment 
of a Morton curve in the following (two examples are depicted in Figure 1). 

We will make use of the following symmetry property of the Morton curve: It 
can be traversed forward or in reverse. The reversal amounts to go through the 
indexing (2) by counting backwards. A quadrant is transformed into the reverse 
ordering by taking the bitwise negation (the one-complement) of its index, 

R{Q) = 2^^ -l-Q. (4) 

Geometrically, this operation mirrors the quadrant around the center point of 
the root cube. 


2.2 The simplicial Morton curve 

The tetrahedral Morton (TM) SFC applies to triangular and tetrahedral red- 
refinement of a mesh (and, conceptually, to higher dimensional simplices) [7]. 
We encounter 1:4 rehnement in 2D and 1:8 refinement in 3D [5], which means 
that the quad-/octree interpretation is still valid. We compute the SFC in a 
bitwise fashion that is an extension of the traditional Morton curve. In order to 
define the TM SFC we introduce the concept of the type of a simplex. 

Definition 3. We describe a d-dimensional simplex T <zW^ by d-\-l ordered ver¬ 
tices xo, ..., Xd G and write T = [xq, ..., x^]. By Xij we denote the midpoint 
between Xi and Xj. 

Bey’s red-rehnement rule [5] for a triangle (d = 2) or tetrahedron (d = 3) 
amounts to dividing the parent simplex T = [xq, ..., x^] into 2'^ subsimplices that 
are defined and enumerated as follows (see also Figure 2): 


— 9 ■ 

To :— [xo,xoi,xo2] 

Ti 

[2:01,2:1, X12], 


T2 := [XQ2,X12,X2] 

T3 

:= [2:01,2:02,2:12], 

To 

:= [xo,xoi,xo2, 

2:03], 

T4 := 

[xoi, X02 

2:03,2:13], 

Ti 

:= [xmiXi^xi2, 

2:13], 

Ts := 

[xoi, X02 

2:12,2:13], 

T 2 

'■= [x02,Xi2,X2, 

2:23], 

Te := 

[2:02,2:03 

2:13,2:23], 

Ts 

■= [2^03,2:^13,3:23 

,2:3], 

T7 : = 

[xo 2 ,2:12 

2:13, 2 : 23 ]- 


(5a) 


(5b) 
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Figure 2: Left: the refinement scheme for triangles in two dimensions. A triangle 
T = [xo,xi,X 2 ] C is refined by dividing each face at the midpoints Xij. We obtain 
four smaller triangles, all similar to T. Right: the situation in three dimensions. If we 
divide the edges of the tetrahedron T = [ccq, xi, 2 : 2 , xa] C in half, we get four smaller 
tetrahedra (similar to T) and one octahedron. By dividing the inner octahedron along 
any of its three diagonals (shown dashed) we finally end up with a partition of T into 
eight smaller tetrahedra, all having the same volume. The refinement rule of Bey is 
obtained by always choosing the diagonal from X 02 to X 13 and numbering the corners of 
the children according to (5). 


Definition 4 (Type of a simplex). We start with a unit square/cube divided as in 
Figure 3 and pick any of the triangles/tetrahedra as root simplex for refinement. 
Each subsimplex in a uniform level L refinement of this root simplex is contained 
in a sub square/subcube of level L and is exactly one of the two (2D) or six (3D) 
simplices from Figure 3. It thus has a unique number, which we define as the 
type of the simplex. 

We start on level 0 with the root simplex T^, which can have any of the 
possible types. In our implementation we pick 0 as the type of the root simplex. 
The TM code m{T) for a descendant T of the root simplex is the interleaving 
of its anchor (lower left) node coordinates with the types of all of T’s ancestor 
simplices [7]. It creates a total order between all simplices of a given level and 
thus establishes the SFC. Here we give a second, recursive definition of the SFC 
that is more suitable for our purposes. 

By Proposition 18 in [7] we obtain one permutation ab G each possible 

simplex type b. It relates the ordering of its children to the SFC such that for 
any d-simplex T with type(r) = b 

m{T„,(o)) < < ■ ■ ■ < m(r^^(2d_i)). (6) 

It places the child Tj in Bey’s order at SFC position ab{i). 

Definition 5. Let T be a level L descendant of such that T's parent P has 
type b and T is the i-th child of P according to Bey’s order (5), 0 < f < 2^. IFe 
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Figure 3: We can separate an axis aligned cube into subsimplices by dividing it along one 
diagonal. We enumerate the resulting subsimplices and call the number of a subsimplex 
its type. Left: We divide a 2D square into two triangles, the lower right one has type 
zero and the upper left one has type 1. Right: We divide a 3D cube into six tetrahedra 
and enumerate them counterclockwise from zero to five (exploded view). 
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Table 1: The local index of a d-simplex T. For each b = type(T), the 2“^ children 
To,... ,T 2 d_i of T can be ordered according to their TM-indices. The position of the 
i-th child according to this order is the local index /loc(Ti). 


call the number crf,(i) the local index of the d-simplex T and use the notation 

hociT) := ab{i). (7) 

By definition, the local index of the root simplex is zero, I\oc{T^) := 0. Table 1 
lists the local indices for each parent type. 

Thus, we know for each type 0 <b < dl how the children of a tetrahedron of 
type b are traversed. This gives us an approach for describing the SFC arising 
from the TM-index in a recursive fashion [12]. By specifying for each possible 
type b the order and types of the children of a type b simplex, we can build up 
the SFC. In Figure 4 we describe the SFC for triangles in this way. In three 
dimensions it is not convenient to draw the six pictures for the different types, 
yet the SFC can be derived similarly from (5) and Table 1. 

Remark 6. In 2D, we will make use of a symmetry property similar to (4).' 
Reversing the TM curve in a uniform refinement of a type 0 triangle results in 
the (forward) TM curve for a type 1 triangle, and vice versa. 
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type = 0 


type = 1 



Figure 4: Left: Using the notation from [12] we recursively describe the space-filling 
curve arising from the TM-index for triangles. The numbers inside the child triangles Ti 
are their local indices /loc(Ti). We write R for the refinement scheme of type 0 triangles 
and F for type 1 triangles. This pattern can be obtained from (5) and Table 1. Right: 
the SFC for a uniform level 3 refinement of the root triangle. 
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Figure 5: We show two L = 2 examples of Morton curve segments that begin with the 
first subquadrant of a tree (left) and end with its last subquadrant (right), respectively. 
In both cases the segment covers one face-connected subdomain. 


3 Illustrated proofs for d < 3 

This section is devoted to proofs that use geometric intuition in two and three 
dimensions. For the cubical Morton curve, the idea is not new (although the 
execution in 3D seems to be). For the tetrahedral Morton curve, this is the 
first such study as far as we know. For abstract proofs for cubes of arbitrary 
dimension d we refer the reader to Section 4. 

3.1 The cubical case 

In this section we prove a set of statements for cubes up to three dimensions 
by providing selected illustrations and covering all possible cases. A similar 
argument has been explored before in two dimensions [3], while an abstract 
proof for two dimensions can be found in [9]. 

We begin with statements that assume a curve that either begins with the 
first subquadrant of the unit cube or ends with its last subquadrant. In a second 
step, we use these statements to prove the final result. All statements are stated 
for arbitrary levels of refinement L > 0. In fact, all statements are trivially true 
for one dimension d = 1 (with no jumps at all); in this section we cover d = 2 
and d = 3. 
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Figure 6: Proof of Proposition 7: The three non-trivial cases that occur in two dimen¬ 
sions (we choose L = 3). The letter F designates a fully covered subtree of level L — 1. 
It is crucial that the lower left corner of the hatched area touches at least one of the full 
subtrees across a face. 



Figure 7: Proof of Proposition 7: Selected cases in three dimensions (j = 1,4 out of the 
seven non-trivial ones). The full subtrees are shaded lightly. Again we exploit the fact 
that the lower left front corner of the last non-empty subtree (hatched) connects to at 
least one full subtree with a lower subtree index across a face. 


Proposition 7. In a quadtree (or oetree) T that is uniformly refined to level L, 
a contiguous segment of a Morton curve that begins with the first subquadrant in 
T creates exactly one subdomain of face-connected quadrants, no matter where 
it ends. 

Corollary 8. In the situation of Proposition 1, a contiguous segment that ends 
with the last subquadrant in T creates exactly one face-connected subdomain, no 
matter where it begins (see Figure 5 for an illustration). 

Proof. Assuming that Proposition 7 is true, we can use the symmetry of the 2 ;- 
curve with respect to reversal to transform the present problem into the setting 
covered in Proposition 7. □ 

Proof of Proposition 7. We proceed by induction over L. Starting with L = 0, 
we only have one element and the statement is true. Supposing L > 0, we can 
identify the number j G [0, 2^) that designates in which level 1 subquadrant 
of the tree the last level L subquadrant of the segment lies. If j =0 then 
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Figure 8: The cases 1-3 and 3-6 in the proof of Proposition 9 for d = 3 dimensions. 
Each of these examples produces two distinct face-connected subdomains. 


the whole segment is contained in a level L — 1 subtree and we can apply the 
induction assumption. Each of the remaining cases produces j full subtrees and 
one possibly incomplete one. That last subtree necessarily contains its first level 
L subquadrant q. Since this subtree produces one subdomain by induction, we 
are done by arguing that the full subtrees are face-connected to each other and 
to q, directly or indirectly. For two dimensions we show the three possible cases 
in Figure 6, all of which satisfy the statement. For three dimensions we proceed 
by enumeration as well; we show selected situations in Figure 7 to conclude the 
proof. □ 

Now that we have identified situations that produce one subdomain only, we 
can prove the main statement for arbitrary segments by a divide-and-conquer 
approach. 

Proposition 9. In a quadtree or octree that is uniformly refined to level L, 
a contiguous segment of the Morton curve creates no more than two distinct 
face-connected subdomains. 

Proof. We proceed by induction over L. Again, the case L = 0 leaves nothing to 
prove. If the segment of the curve is contained in one level L—1 subtree, the proof 
is finished by induction. Else we have one subtree in which the segment begins, 
zero to 2^^ — 2 fully covered subtrees, and one subtree in which the segment ends. 
To the first nonempty subtree we can apply Corollary 8, while Proposition 7 
applies to the last one. Thus we know that the possibly incomplete subtrees 
lead to one connected piece each. The case of two nonempty subtrees is thus 
completed and it remains to consider three or more. 

Now, whenever any two adjacent nonempty subtrees have even-odd numbers, 
they are face-connected since at least one of them must be full. This covers the 
remaining three- and four-subtree cases in two dimensions. In three dimensions, 
this clears all situations with three non-empty subtrees. Since we can further 
reduce the number of remaining cases by symmetry, it remains to examine the 
subtree ranges {i,... ,i + 3) through (i,..., 7) for i = 0,...,3. All of these cases 
satisfy our claim; we illustrate a few in Figure 8. □ 
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Figure 9: Left: a segment of the 2D SFC on a level 4 refinement of T 2 with six face- 
connected components (shaded pink). The number of face-connected components in 2D 
can be as high as 2{L — 1); this estimate is sharp. Right: a 3D level 2 refinement of Tg 
with four (= 2L) face-connected components. We prove that an upper bound on the 
number of face-connected components is 2L + 1 and conjecture that 2L is sharp. 


We have completed the necessary proofs for a uniform space division into 
cubes in d < 3. In Section 4.1 we extend the proof to arbitrary dimension d. 
The case of adaptive space divisions is considered in Section 4.3. 

3.2 The simplicial case 

In this subsection we examine the number of face-connected components of a 
segment of the tetrahedral Morton SFC, d = 2 or 3. As we show in Figure 9, 
there exist cases where the number of face-connected components in a uniform 
2D level L refinement can be as high as 2{L — 1). We show that this is in fact 
a sharp upper bound. We also show that in three dimensions the number of 
face-connected components does not exceed 2L -|- 1. There exists an example 
with 2L face-connected components and we conjecture that 2L is in fact the 
sharp estimate. The proof of these bounds is fairly analogous to the results for 
cubes and relies and a divide-and-conquer approach by splitting the segment 
into subsegments of which we know the number of face-connected components. 
The main difference to the cubical Morton curve is that we do not have a strong 
symmetry property like (4), and thus an analogue to Corollary 8 only exists in 
a weaker form. 

Lemma 10. The following two properties hold for the TM-index in 2D, where 
we consider a uniform level L refinement of an initial type 0 triangle T. 

• Each type 1 subsimplex is face-connected to a type 0 subsimplex with a 
greater TM-index. 

• Each type 0 subsimplex that is also a descendant of the level 1, type 1 
subtriangle T 3 is face-connected to a type 1 subsimplex with a greater TM- 
index. 

Proof. The respective face-neighbor is the top face-neighbor for the type 1 sub¬ 
simplex and the face-neighbor along the diagonal face for the type 0 subsimplex; 
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Figure 10: Illustration of Lemma 10. In 2D, choose any subsimplex 5*. If its neighbor 
along the top face S'' is inside the root triangle, then m(S*) < m(S'). This condition is 
always fulfilled by any type 1 triangle and by type 0 triangles that are descendants of 
the middle level 1 subtriangle. 


see Figure 10. For type 0 we additionally require that the subsimplex is a de¬ 
scendant of Ts, since this ensures that the face-neighbor along the diagonal face 
is inside the root triangle. Despite this detail, the proofs for both items are 
identical, and we only present one for the first. 

Let S denote an arbitrary type 1 subsimplex of level L and let S' be its 
neighbor across the top face. If S and S' share the same parent P then there 
are two cases, which we also see in Figure 4: Either type(P) = 0, then the local 
index of S' is 2 and that of S' is 3, or type(P) = 1, in which case the local index 
of S is 0 and that of S' is 1. Thus, in both cases the TM-index of S must be 
smaller than that of S'. We suppose now that S and S' have different parents, 
which implies L > 2, and denote these different level L — 1 subsimplices by P 
and P'. The only possible combination is that type(P) = 1 and type(P') = 0, 
and that P and P' are neighbors along P’s top face. Therefore, by an induction 
argument, m{P) < m{P'), and since the TM-index preserves the local order 
under refinement, each child of P has a smaller TM-index than each child of P'. 
In particular we find m(S) < m{S'). □ 

Let us now show a 2D analogue to Proposition 7. 

Lemma 11. Consider a triangle T that is uniformly refined to level L. If T 
has type 0, then a contiguous segment of the SFC ending in the last level L 
subsimplex has just one face-connected component. If T has type 1, then this 
holds for segments starting in the first level L subsimplex. 

Proof. We only show the statement for type(r) = 0, since we can then use 
the symmetry of the 2D curve (Remark 6) to obtain the result for the case 
type(r) = 1. We proceed by induction over L. 

For L = 0 there is only one possible segment and it is connected. For L > 0, 
let j G {0,1, 2, 3} be the local index of the level 1 subtree T' of T in which the 
first level L subsimplex of the segment lies. If j G {0,1,3}, then the type of T' 
is 0 and the statement follows by induction with the same argument as in the 
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proof of Proposition 7. Thus, let j = 2, i.e., the segment starts in the type 1 
subtree of T. The part of the segment that is not inside T' is the full last subtree 
of T (local index 3) and thus it is face-connected. With Lemma 10 we conclude 
that each subsimplex in the subsegment in T' is face-connected to a simplex 
with greater TM-index. Since this holds particularly for the last subsimplex in 
T', the subsegment is also face-connected to the subsegment making up the last 
subtree. Thus, the whole segment is face-connected. □ 

For all other segments beginning with the first or ending in the last level L 
subsimplex, and notably for all of those segments in 3D, we obtain an upper 
bound of L -|- 1 face-connected components, which we show in the next two 
lemmas. 

Lemma 12. Let a segment of the space-filling curve for a uniform level L refined 
d-simplex consist of several full level 1 subsimplices plus one single level L simplex 
at the end or at the beginning, then this segment has at most two face-connected 
components. 

Proof. Similarly to the last paragraph in the proof of Proposition 9, and in 
analogy to Figure 8, we can show this claim by enumerating all possible cases. □ 

Lemma 13. If a d-simplex is uniformly refined to level L, then any segment of 
the space-filling curve ending in the last subsimplex or starting in the first has 
at most L 1 face-connected components. 

Proof. Consider the case that the segment starts in the first simplex. For L = 0 
there is only one possible segment consisting of the unique level 0 subsimplex 
and it is thus connected. Let now L > 0. Since the segment begins at the very 
first level L subsimplex, we can separate it into two parts. The first part at the 
beginning consists of 0 to 2*^ — 1 full level 1 subtrees, and the second part is one 
possibly incomplete level 1 subtree. 

By the induction assumption, the second part has at most L face-connected 
components. From Lemma 12 we obtain that the first part together with the 
first level L subsimplex of the second part has at most two face-connected com¬ 
ponents. Since this first level L subsimplex is contained in one of the components 
of the second part, we obtain 


L + 2-l = L + l 


( 8 ) 


components in total. 

If the segments ends in the last simplex, the order of parts is reversed. The 
first part of the segment is the part in the level 1 subtree where the segment 
starts, and the second part consists of the remaining full level 1 subtrees. We 
obtain the bound on the number of face-connected components using the same 
inductive reasoning as above. □ 

We have so far argued the connectivity of specific kinds of SFC segments. 
This suffices to proceed to arbitrary segments of the tetrahedral Morton SFC. 
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Proposition 14. Any contiguous segment of the space-filling curve of a uniform 
level L > 2 refinement of a type 0 simplex has at most 2{L — 1) face-connected 
components in 2D and 2L + 1 face-connected components in 3D. For L = there 
are at most two face-connected components (this applies to both 2D and 3D). 

Proof. We first show that for d < 3 the number of face-connected components 
is bounded by 2L -|- 1: If a given segment is contained in a level 1 subtree, we 
are done by induction. Otherwise we can divide the segment into three (possibly 
empty) pieces: First, the segment in one incomplete level 1 subtree ending at 
its last level L subsimplex, then one contiguous segment of full level 1 subtrees 
and finally a segment in one (possibly incomplete) level 1 subtree that starts at 
its hrst level L subsimplex. Lemma 13 implies that the first and the last piece 
have at most L face-connected components each. By Lemma 12, the second 
piece has one or two face-connected components, and if the number is two, then 
it is face-connected to the first or to the third piece. Thus, it adds only one 
face-connected component to the total number, and we obtain at most 

L + 1 + L = 2L + 1 (9) 

face-connected components. 

Let us now specialize to 2D. We conclude from Lemma 11 that the first sub- 
segment only adds more than one face-connected component if it is contained 
in the only level 1 subtree of type 1 (local index 2). Similarly, the second sub- 
segment only adds more than one face-connected component if it is contained 
in a level 1 subtree of type 0. In particular, if both subsegments add more than 
one connected component, the second subsegment is contained in the last level 
1 subtree (local index 2). Thus, the middle subsegment is empty in this case. 

If both of these subsegments have less than L face-connected components, 
there is nothing left to show since the overall number of components is then 
less than or equal to 2 (L — 1). So suppose that one of the subsegments has L 
face-connected components and the other one has at least L — 1. We depict this 
situation in Figure 11. We observe that the first and second level L simplex 
in this first segment are face-connected to the first and second level L simplex 
in the second segment. If, however, the second subsegment has L connected 
components then its last two level L simplices are face-connected to the last two 
level L simplices of the first subsegment. 

We thus can subtract two connected components from the total count, which 
leads to at most 

L + L-2 = 2(L-l) (10) 

face-connected components in total. □ 

We briefly discuss whether we can sharpen these bounds. In 2D, this is not 
possible by counterexample; see Figure 9. In 3D, we construct a segment with 
2L face-connected components using the SFC-indices 22-25 of a uniform level 2 
refinement of a type 0 tetrahedron. We believe that the case that the first and 
the last piece described in the proof of Proposition 14 have L face-connected 
components each and that additionally the middle piece adds one component 
does not occur. 
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Figure 11: An illustration of the 2D case in the proof of Proposition 14 for L = 4. 
The bottom segment has the maximal number of L face-connected components. Since 
its first and second triangle (on the left, shaded in pink) are connected with the top 
segment, the possible number of connected components is reduced by two. If however, 
the second segment has L face-connected components then its last two triangle (on the 
right) are connected with the bottom segment. Thus, the number of face-connected 
components is less than or equal to 2L — 2. 


Conjecture 15. In 3D, the number of face-connected components is bounded by 
2L. This estimate is sharp. 

4 Proofs for arbitrary dimension 

The construction of the cubical Morton curve generalizes readily to arbitrary 
space dimension d. We should suppose that the main result (that any contiguous 
segment consists of at most two face-connected subdomains) generalizes as well. 
Indeed, we propose two different ways to prove this in the following Sections 4.1 
and 4.2. The first is closer to the geometric approach we have been using in 
Section 3, while the second is more formal and paves the way for quantitative 
studies of the frequency of disconnections in Section 5 below. 

We close this section with the extension of the proofs for both cubical and 
tetrahedral curves from uniform to adaptive meshes (see Section 4.3), which is 
the remaining step to establish Theorems 1 and 2, and discuss implications for 
a forest of octrees (Section 4.4). 

4.1 Induction proofs for d-cubes 

We use induction over both the dimension and the level of subdivision to prove 
the main statement for all dimensions d > 0. These proofs imply the statements 
of the previous Section 3.1 as special cases. For convenience we denote any d-tant 
as a quadrant. We make use of the following definition of subtree ranges. 

Definition 16. Let the spaee dimension be d > 0. For any 0 < d' < d and 
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0 < A: < 2'^ , we define the following interval containing 2^' integers, 

if =2'^'[k,k + l). (11) 

We use this interval to denote a specific contiguous range of subtree indiees. 

We define the following auxiliary statements, first considering a one-sided 
segment and then a general two-sided one. 

Proposition 17. If a segment of a Morton curve is fully contained in the level 
1 subtrees enumerated by a given if and contains the first or last subquadrant 
in this range of subtrees, then it corresponds to one face-eonnected subvolume. 

Proof. By symmetry of the Morton curve, we can restrict the discussion to the 
case of the first subquadrant. Let us begin by proving the statement for subdi¬ 
vision level L = 1. By (11) the lowest subtree index in the segment is k2'^ . This 
number has d' zero bits from the right. All other indices in I^ have one or more 
ones in the lower d' bits while being bitwise identical in the higher bits. For any 
of these indices we can flip the low bits to zero one by one, effectively transi¬ 
tioning through face neighbors and monotonously decreasing the index until we 
reach . This whole sequence of face-connected subtrees is contained in I^ . 
In conclusion, all trees in if are face-connected to k2'^' and thus to each other. 

Now let L > 1 and assume the above statement for L — 1. To prove it for 
L we make an induction over d'. If d' = 0 we have a single subtree and can 
readily invoke the induction assumption for L — 1. Else there are two possible 
cases: Either the segment is fully contained in one of or 

apply the induction over d'. Otherwise contains full subtrees only and 

the segment reaches into - Each nonempty subtree j in the latter interval 
must contain its first subquadrant, which has a face connection to the full tree 
j — 2^' G Since by the proof for L = 1 all subtrees in face- 

connected, we are done. □ 

Proposition 18. If a segment of a Morton curve is contained in the level 1 
subtrees if, it produees no more than two distinct face-connected subvolumes. 

Proof. Again let us prove the statement first for L = 1. If d' = 0 we have 
just one level 1 subquadrant that clearly satisfies our claim. For positive d' we 
distinguish the following cases. If the segment is fully contained in either I^jT^ 
or / 2 fe+i’ s-PPly the induction on d'. Else we know that the last subquadrant 
of l 2 ]fi^ the first of ™ segment. By Proposition 17 we have at 

most two disconnected pieces and the statement holds. 

If L > 1 the case d' = 0 reduces to the same statement for L — 1 and we are 
done by applying the induction over L. Else, the proof proceeds unchanged as 
above with the desired result. □ 

We have implicitly proved the main result for any uniform level L subdivision, 
since a level 0 subtree trivially satishes our claim, and otherwise the root cube 
is the union of the level 1 subtrees if 
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4.2 A non-inductive proof for d-cubes 

In this section we elaborate on the formalism of the Morton index (see Sec¬ 
tion 2.1) to obtain the result without induction. The tool we use is the map 11 
from the index Q = {q^... q ^)2 to a subset of stated in (3). For 1 < r < d, 
we define the coordinate along axis r, 

Qr = {ql---qr) 2 - (12) 

The map 11 (Q) may be written as 

11(Q) = [2-^Qi,2-^(Qi + 1)] X • • • X [2-^Qd, 2-^{Qd + 1)]. (13) 

Lemma 19. If Q and Q are such that Qk < Qk for all 1 < k < d, then Q < Q. 
Proof. The order of the bits in Qk (Qk) is the same as their order in Q (Q). □ 

Theorem 20. For any index Q , the interior o/Tq = Uq=o^(^) star-shaped 
(and thus Yq is face-connected and contractible). 

Proof. If (5 G {0,..., then so are all Q such that Qfc ^ Qk for all 1 < 

k < d. The domains of these quadrants define a box between the origin and the 
corner of ^{Q) farthest from the origin, 

B{Q) := [0, 2-^(Qi + 1)] x • • • x [0, 2-^(Qd + !)]• (14) 

Indeed, Tq is the union of these boxes, Yq = Ug^Qi?(Q), and the union of their 
interiors is the interior of Yq. Each of these boxes contains the midpoint of H(0) 
in its interior and is star-shaped with respect to it. Therefore the interior of Yq 
is star-shaped with respect to that point as well. □ 

Corollary 21. = Q (\ means bitwise-or) for all Q G ..., 

then the interior ofY = DQ^Q^tart^yQ) is star-shaped. 

Proof. If = Q, then the 1-bits of Q—Q®*®'*'* are a subset of the 1-bits of Q. 

One can then verify that {Q — = Qk — for all 1 < A; < d. Therefore 

U JjQstanll(Q-Q"*“*) is Y translated by the vector ..., -2-^Qf^^^). 

^end_^start 

This is the same as ^q=q ^ ^{Q}^ which is star-shaped by Theorem 20. □ 

Corollary 22. If = Q (k. denotes the bitwise and operator) for all 

Q G {(5®*“'’*,..., then the interior ofY is star-shaped. 

Proof. Mirroring every quadrant about the midpoint of the unit cube does not 
change the shape of Y. The mirror of 11(Q) is 11(12(0), where R{Q) denotes 
the bitwise negation (4). Therefore R{Q)\R{Q'^'^‘^) = i2(Q&0°‘^) = R{Q) for ah 
R{Q) G {12(Q'^°‘^),..., 12(Q"*”*)}. □ 

Theorem 23. The interior of Y is star-shaped, or Y is the union of two sets 
whose interiors are star-shaped. 

Proof. Let q be the most significant bits common to all of ..., Q™^}. We 

can split the segment into > (^011... 1 ) 2 } and {(glOO ... 0 ) 2 ,. • •, 0“^}- 

The interior of the domain of the first segment is star-shaped by Corollary 22; the 
interior of the domain of the second segment is star-shaped by Corollary 21. □ 
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4.3 From uniform to adaptive meshes 

We have completed the necessary proofs for a uniform space division, in the case 
of cubical refinement for any space dimension d, and previously for triangular 
and tetrahedral refinement (see Section 3.2). As we state in this section, an 
adaptive space division does not require any more effort (see also [3, page 176]). 

Proof of Theorems 1 and 2. Any adaptive tree of quadrants with level < L can 
be refined into level L quadrants exclusively. This operation does not change the 
connectivity between boundaries of the designated subdomain. In particular, the 
number of face-connected subdomains remains unchanged and the proof reduces 
to applying Propositions 9, 14 (only d < 3) or 18 (any d) above. □ 

4.4 From one tree to a forest 

If we consider a forest of octrees as in [4,8,19], a contiguous segment of the 
Morton curve may traverse more than one tree. In this case, the segment nec¬ 
essarily contains the last subquadrant of any predecessor tree, as well as the 
first subquadrant of any successor tree in the segment. In the cubical case, we 
know by Proposition 17 that no jumps can occur at all (when not counting the 
transition between two successive trees as a jump). For the simplicial case, we 
may use Lemmas 11 and 13 to obtain the bounds L -|- 1 (2D) and 2L -|- 1 (3D). 

5 Enumeration of face-connected segments 

We would like to examine not only how many pieces an SFC segment can have, 
but also how frequently segments of different numbers of pieces occur. To this 
end, we propose a theoretical lower bound for the cubical case and supply nu¬ 
merical studies for both cubical and tetrahedral SFCs. 

5.1 Lower bound on fraction of continuous segments 

Theorem 24. The fraction of continuous segments of length I of the level-L 
d-dimensional, eubical Morton eurve is 

4>d,L,i > 2d _i - 

Proof. Let s = ..., — 1} be the first discontinuous 

segment of length 1. As in the proof of Theorem 23, s divides into two continuous 
segments, {Q®*®'*'*,..., (qOll... 1 ) 2 } and {(glOO ... 0 ) 2 ,..., where q are 

significant bits that are common to all numbers in the segment. By the same 
reasoning as used in Corollary 21, the shapes of the two segments are not affected 
by q: changing q translates the whole domain. Therefore, as this is the first 
discontinuous segment, q must be (0... 0), i.e., the two continuous pieces of the 
segment are s~ = ..., 2*^ — 1} and = {2 ^,..., for some k. Since 

cannot have a more significant 1-bit than 2^, < 2^+^ — 1. We note 

that bitwise negation of the first A: -|-1 bits, R{Q) '■= — 1 — Q, induces a map 
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^{Q) ^{R{Q)) that is the reflection about the midpoint of the box formed 

by UQ^(i"^o(g). 

We first want to find a lower bound for We note that 2^ G s"*", so if 

h2(Q*) is a face neighbor of 0(2^), then Q* 0 s“, and if Q* <2^, then because 
of the definition of s~, Q* is less than every Q ^ s~, including Let 

j = k mod d and let j' = d — j — 1. The Q* that plays the crucial role is the 
neighboring quadrant that is closer to the origin in the {j + l)th direction. We 
have to subtract one from the {j + l)th coordinate of 2^: i.e., (2^)j-|-i = in 

the notation of (12), so Qj^i = — 1, while the other coordinates are the 

same, viz. zero. Therefore 


[fc/dj-times 

Q * = (0 ... 0 ( oT ~ o 1 ' o ]^)2 

/-times /times 

[fc/dj 


2 = 1 


2/^2 


—id 


= 2'' 


S 2 

2=1 2=[A:/dJ-|-l 

- oo 

-2^' ^ 2"*'^ 


—id 


2 ^ - 1 


2 = 1 


[k/d] 

2=1 


[k/d\ 


2^Y,2 


—id 


= 2 ' 


2 = 1 

^ OO 

—2^ y 2 

2d _ 1 

i= [k/d] H-l 


(16) 


—id 


2^ - 2 ^ 
2'^ - 1 ■ 


By the dehnition of there is a continuous segment of length I that 

starts at each Q G {0,..., Q*}. For each of these continuous segments, there is 
another continuous segment, obtained by the reflection map R, than ends with 
Q G {i?(Q*),..., 2^~^^ — 1}. We want to show that these two sets of continuous 
segments are distinct, i.e., that there is no segment of length I that starts with 
g < g* and ends with Q + I — 1 > R{Q*). We thus have to show that the 
shortest segment with endpoints in each set, {g*,... ,i?(g*)} is longer than I, 
i.e., I < R{Q*) -Q* + 1 = 2^+1 - 2gL 

We will prove this bound by finding an upper bound for g®“'^. We note that 
R(2^) = 2^ — 1 £ s~, so fl{R{Q*)) is a face neighbor of 12(2^ — 1), so by the 
same reasoning as above, R{Q*) must be greater than and thus 

I = g«>^d - g®t^rt + 1 < {R{Q*) - 1) - (g* + 1) + 1 = 2^+1 - 2Q* - 2. (17) 


We have shown that there are at least 2(g* + 1) continuous segments in the 
first 2^^^ segments, each of which begins and ends in the range {0,..., 2^^^ — 1}. 
None of the numbers in these segments has more than A; + 1 significant bits, so 
by the same reasoning as in Corollary 21, adding a multiple of 2^~^^ to each 
number in one of these segments is a translation of its domain, and is thus 
continuous. Therefore there are at least 2(g* + 1) continuous segments for every 
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2^+1 segments, and thus 


4‘d,L,l ^ 2 ^“*"^ 2 ^ 


2{Q* + 1) _ Q* + 1 


2fc(^2“-lj 2 ^ 

2k _ 2i + 2'^ - 1 
2^(2'=* - 1) 



(18) 


- 2k{2^ - 1 ) 


{j < d hy definition) 


1 


2'^ - 1' 


□ 


5.2 Computational studies—cubical Morton curve 

Having shown that a segment of a Morton curve is composed of one or two 
face-connected subdomains, a natural question to ask is how many of each type 
there are. More formally, we ask: for a given dimension d, recursive level L, and 
segment length I, what fraction (jid^L.i of the 2'^^ — (/ — 1) possible segments are 
in one face-connected piece? 

This question can be answered recursively. Each segment of length 1 > 1 on 
level L can refine to 2^^ segments on level (T-l-1) with lengths between 2'^(/—2)-|-2 
and 2 '^l, as illustrated in Figure 12. We divide connected segments into two 
categories: weakly connected, when the first and last quadrants in the segment 
are (face-)adjacent, and strongly connected, when they are not. Disconnected 
segments only refine to disconnected segments. Strongly connected segments 
only refine to strongly connected segments. Weakly connected segments refine 
to all three types: how many of each depends on the direction in which the first 
and last quadrants are adjacent. 

We give pseudocode for this recursive calculation in the function Enumerate 
(Algorithm 1). This algorithm is implemented in the Python script morton.py.^ 

Enumerate calls on some lookup tables: RefineOne(d,l) (Figure 13) counts 
how many disconnected, strongly and weakly connected segments of length I 
are refined from one d-dimensional quadrant (the weakly connected segments 
are broken down by the direction in which the first and last quadrant are ad¬ 
jacent); RefineWeak(d,j,r) (Figure 14) counts how many disconnected, strongly 
and weakly connected segments with r quadrants in the end-families are refined 
from one weakly connected segment in direction j (a weakly connected segment 
only produces weakly connected segments in the same direction). 

In Figure 15, we use Enumerate to calculate the fraction of connected seg¬ 
ments 4>d,L,i for d = 2 and d = 3 for large values of L. We observe that (j)d,L,i 
tends to vary between 1/2 and 1/(2'^ — 1). (We had proved in Section 5.1 that 
the latter is indeed a lower bound.) 

^https: //github. com/cburstedde/p4est/tree/develop/doc/morton/mortoii.py 
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Figure 12: The refinement of disconnected (left), weakly connected (middle), and 
strongly connected segments (right). Each coarse segment (top) refines to one of 2^'^ = 16 
possible refined segments (bottom), each starting with s G {sq, •.., S 3 } and ending with 
e G {eo,... , 63 }. Disconnected segments refine to disconnected segments; strongly con¬ 
nected segments refine to strongly connected segments; weakly connected segments refine 
to disconnected (e.g., { 53 ,..., eo}), weakly connected (e.g., {s 2 ,.. •, eo}), and strongly 
connected segments (e.g., {sq, ■ ■ ■, es}). 


r ■ ■ ■ ■ 

52 ^; 


L . . . . 

““51 


RefineOne(2,l) 

Us = 4: ({so}, {si}, {S2}, {S3}) 

RefineOne(2,2) 

nd = l ({si, S2}), n^,i = 2 ({so, si}, {s2, S3}) 

RefineOne(2,3) 

nw ,2 = 2 ({so, Si, S2}, {Si, S2, S3}) 

RefineOne(2,4) 

Us = 1 ({so,Sl,S2,S3}) 


Figure 13: We list the RefineOne(fi,/) tables used in Enumerate (Algorithm 1) for d = 2 
as an example. Unlisted values are zero. 


Refine Weak (2,2,2) 

= 1 ({S3, ■ • •, 

eo}) 



Refine Weak (2,2,3) 

^iu,2 — 2 ({ 52 , • • 

■, eo}, {s3, • 

■•,ei}) 


Refine Weak (2,2,4) 

Us = 3 ({si,..., 

eo}, {s 2 , • • ■ 

, ei}, {s 3 ,. 

..,e2}) 

Refine Weak (2,2,5 < r < 8) 

ns = 9-r ({so, 

• • ■ , Cr-s},- • 

• ,{S8—r, • ■ 

,e3}) 


Figure 14: We list the RefineWeak((i,j,r) tables used in Enumerate (Algorithm 1) for 
d = 2 and j = 2 as an example. The start points and end points refer to Figure 12 
(middle). Unlisted values are zero. 
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Algorithm 1: Enumerate (d, L, 1) 


Data: dimension d > 1, level L > 0, segment length 1 < I < 2^^. 

Result: (n^, rig, n^^i, ..., n^^d)i the number of segments of length I that are 
disconnected, strongly connected, and weakly connected in each 
direction. 

1 if Z = 1 or L = 0 then [define segments of length 1 to be strongly connected] 

2 I return (0, 2^^, 0 ,..., 0) [one segment for each quadrant] 

3 end 

4 {jld^ '^s-i • • • 1 '^w,d') ^ (0; • ■ • ; 0) 

5 c \l/2‘^~\ [compute the shortest length that can refine to length 1] 

6 if I mod 2‘^ <2 then C •(— 1 + [l/2‘^\ 

7 else 2 + [Z/2^J 

[... longest ... ] 


8 for k G {c,..., C} do 

9 if fc = 1 then 


10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 
23 


N G- 2'^^^ of coarse quadrants] 

{ud, Us, riu,,!,..., riyj^d) -G- {ud, Us, • ■ ■, + N* RefineOne (d,/) 

else 

r ^ I — {k — 2)2‘^ yf of children in two coarse end quadrants] 

m G- min{(r — 1), 2'^+^ — (r — 1)} [# of ways to split between two families] 
(Arf, Ns, iVu,,i,..., N,jj^d) G- Enumerate (d, L - 1, fc) 
nd <— nd + Nd * TO [disconnected -G disconnected] 

Us <— rig + Ns * m [strongly connected -G strongly connected] 

for 1 < j < d do 

[ (ud, Ug, riujj) t- (ud, Ug, n^, j) + N^^j* RefineWeak (d, j, r) 

end 

end 

end 

return (rid, rig, ..., n^^d) 



Figure 15: We plot the fraction of continuous segments of length I, (j)d,L,i, for d = 2 and 
L = 30, (left) and d = 3 and L = 21 (right), for one thousand log-uniformly randomly 
sampled lengths. 
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Figure 16: The relative count of SFC segments by number of connected components 
and the average length (right y-axis) of these segments. We exclude the segments of 
length 1. Left: the distribution for a uniform level 8 refined triangle. We observe 
that almost 98% of all SFC segments have three connected components or less. 63.7% 
are connected, 29.7% have two connected components and 4.4% have three connected 
components. Right: the distribution for a uniform level 5 refined tetrahedron. Here, 
more than 93% of the segments have three connected components or less with 61.0% 
having exactly one connected component, 22.1% with two connected components and 
10.7% with three connected components. The highest number of segments occuring 
agrees with Proposition 14 (2D) and Conjecture 15 (3D). 


5.3 Computational studies—simplicial Morton curve 

For simplices, we enumerate all possible SFC segments for a given uniform re¬ 
finement level and compute the number of their face-connected components. We 
achieve this by performing a depth-first search on the connectivity graph of the 
submesh generated by the segment^. 

While we investigate the fraction of connected SFC segments of a given par¬ 
ticular length I among all segments of length I in Section 5.2, here we we compute 
the fraction of connected segments of any length among all possible segments. 
More precisely, we compute for each possible count of connected components 
the chance that any randomly choosen SFC segment (with a random length) has 
exactly this number of connected components. 

For a uniform level 5 refined tetrahedron we obtain that 61% of all SFC seg¬ 
ments are connected and only 7% do have four or more connected components. 
For a uniform level 8 refined triangle, about 64% of the segments are connected 
with 2% of the segments having four or more components. For cubes and quadri¬ 
laterals, the respective ratios of connected segments are 60% and 71% (here we 
know that the disconnected segments have exactly two components). We collect 
these results in Figure 16 and Table 2. 

^https://github.com/holke/sfc_conncomp 
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Level 5 

Level 8 


Quads 

Cubes 

Triangles 

Tets 

Quads 

Triangles 

Connected 

Non-connected 

71.5% 

28.5% 

60.0% 

40.0% 

63.8% 

36.2% 

61.0% 

39.0% 

71.4% 

28.6% 

63.7% 

36.3% 


Table 2: The relative counts of connected and non-connected segments across all possible 
SFC segments of a uniform level 5 and level 8 (2D only) refinement, excluding segments 
of length one. Here we average over a uniform distribution of lengths (in contrast to the 
log-uniform distribution used in Figure 15). 


6 Conclusion 

We prove in this document that the classical Morton or z-curve does not lead to 
a fragmentation of the root cube into more than two disconnected subdomains. 
Its loss of continuity in comparison to the Hilbert curve is thus controlled. This 
is in line with experimental results that establish the suitability of the Morton 
curve for numerical applications. 

We show that the bound for the recently proposed tetrahedral Morton cube is 
of order L and thus growing with the level of refinement. Yet, we can demonstrate 
numerically that the fraction of connected to non-connected segments is close to 
the cubical case. In practice, we may expect both approaches to behave similarly. 

Our result would appear relevant to make informed choices about the type 
of space filling curve to use, for example in writing a new element-based parallel 
code for the numerical solution of partial differential equations, or any other code 
that benefits from a recursive subdivision of space. Our theory and experiments 
support the existing numerical evidence that a fragmentation of the parallel 
partition is not observed. 
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